TO ALL WHOM IT MAY CONCERN: 



BE IT KNOWN THAT WE, TAKAYUKI NISHIMURA, a 
citizen of Japan residing at Tottori, Japan, YASUYUKI 
NOMIZU, a citizen of Japan residing at Kanagawa, Japan, 
HIROYUKI SAKUYAMA, a citizen of Japan residing at Tokyo, 
Japan, JUNICHI HARA, a citizen of Japan residing at 
Kanagawa, Japan, NEKKA MATSUURA, a citizen of Japan 
residing at Kanagawa, Japan, TAKANORI YANO, a citizen 
of Japan residing at Kanagawa, Japan, TAKU KODAMA, a 
citizen of Japan residing at Kanagawa, Japan, TOSHIO 
MIYA2AWA, a citizen of Japan residing at Kanagawa, 
Japan and YASUYUKI SHINKAI, a citizen of Japan residing 
at Kanagawa, Japan have invented certain new and useful 
improvements in 

METHOD FOR REVERSIBLY TRANSFORMING DATA FORMAT, IMAGE 
PROCESSING APPARATUS, PROGRAM FOR REVERSIBLY 
TRANSFORMING DATA FORMAT, AND COMPUTER-READABLE 
RECORDING MEDIUM THEREOF 



which the following is a specification:- 
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BACKGROUND OF THE INVENTION 

1. Field of the Invention 

The present invention generally relates to a 
transforming process for transforming a format of image 
data, and in particular, a method for reveraibly 
transforming data format, an image processing apparatus, 
a program for reversibly transforming data format, and a 
computer-readable recording medium thereof, in which a 
color space transforming process is conducted in a system 
where a lossless/lossy compression/extraction is conducted 
with respect to the image data in a color space. 

2. Description of the Related Art 
A color space is one of areas in a vector space 
that is a three dimensions or higher dimensions . For 
example , a color coordinate system is defined by three 
dimensional vectors being three linear independences as 
bases. In the color coordinated that is generally used, 
each of R (Red) , G (Green) , and B (Black) are defined by 
a central wavelength . When a certain three dimensional 
color coordinate system is given, any three dimensional 
linear color coordinate other than the certain three 
dimensional color coordinate can be expressed by a 
reversible 3x3 matrix (non-singular) , 

In practice, various color coordinates (color 
spaces) are defined by various reasons. For example, in 



3 



a case of displaying data at a monitor, a RGB coordinate 
system is used in almost all digital color images. For 
example, a fixed range being 8 bits /coordinate is 
preferable. In a case of a compression requiring a color 
noncorrelativity, the RGB coordinate system is not suitable 
but another color coordinate system such as a YIQ coordinate 
system is more suitable. Also, a YUV coordinate system and 
a YCrCb coordinate system are suitable. Regarding a 
separation between a brightness and a chrominance relating 
to a luminosity, these opponent color coordinates systems 
attempt to improve a separation of a characteristic and an 
appearance of a change of the chrominance in the s ame 
brightness 

In addition, in image for a print, a subtractive 
color system such as a OCT (Cyan, Magenta, and Yellow) is 
used. Alternatively, a perfect four dimensional space such 
as a CMYK (Cyan, Magenta, Yellow, and Black) is used. 

The data compression is extremely useful to 
store and transmit a large amount of data. For example, 
a necessary time to transmit image data such as a facsimile 
transmission of a document can be significantly reduced by 
compressing the image data and reducing the number of bits 
necessary to reproduce an image. 

Conventionally, there are various data 
compressing methods . The compressing methods can be 
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classified into two categories : that is, a lossy coding and 
a lossless coding. In the lossy coding, loss of image data 
maybe occur. Accordingly, an original image can not be 
guaranteed to be produced perfectly. An purpose of the 
lossy compression is to prevent a change caused by 
reproducing the image data from being unpleasant and 
conspicuous. In the lossless compression, all the image 
data are stored and compressed in a state in that the image 
data can be perfectly 'reproduced. 

In the lossless compression, an input symbol 
or brightness data are converted into an output code word. 
As input data, image data, voice data, one dimensional data 
(for example, data changing temporally) , two dimensional 
data (for example, data changing to two space axis 
directions) , and multi -dimensional /multi- spectrum data 
are used. When data are successfully compressed, the code 
words of the compressed data are expressed by the n umb er 
of bits fewer than the number of bits of the input symbol 
(or the brightness data) . As the lossless coding method, 
a dictionary coding method (for example, Lempel-Ziv coding 
method) , a run length coding method, an arithmetic coding 
method, and an entropy coding method are used. In a case 
of the lossless image compression, the basic of the 
compression is a prediction, a context, and a coding, a 
JBIG standard for a facsimile compression and a DS>CM 
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(Difference Pulse Coding Modulation that is an option of 
a JPEG standard) for a continuous -tone image are examples 
of the lossless compression for the image. In the lossy 
compression, the input symbol and the brightness data are 
converted into output code words after quantized. Purposes 
of the quantization is to store characteristic data being 
important portions of data and also to eliminate 
unimportant data. In the lossy compression, a conversion 
for an energy concentration is utilized before the 
quantization. Abase line JPEG is one example of the lossy 
coding method for the image data. 

Conventionally , a color coordinate 
transformation has been utilized for the lossless 
compression with the quantization. In the lossless 
compressing system or the lossless/lossy compressing 
system, a main essential requirement is reversibility and 
effectiveness of a transformation. In other 
lossless/lossy compressing systems, a non-correlation of 
colors is one factor other than the effectiveness of a 
reversible transformation- For example, the 3x3 matrix is 
used for the lossy compression only. Because factors are 
not integers, when the non-correlation is required, errors 
are involved in a process while data is alternately 
compressed and extended. 

When the color space is transformed, a 
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calculation accuracy becomes a problem* For example, when 
8 bits are input, 10 or 11 bits are required to transform 
the color space. Higher accuracy is required for an 
internal calculation, so that a stable color space can be 
obtained . Errors are accumulated in a result after an image 
is converted and compressed from the RGB color space and 
then extended to restore the image in the RGB color space 
even if a sufficient accuracy is applied to processes being 
repeated. Accordingly, a final color does not match the 
same as an original color. This is called color shift that 
occurs since the color space is unstable, 

Methods are provided to realise a reversible 
color space transformation by an integer operation having 
predictable accuracy (for example, refer to IDS or 
Cross-References 1 and 2) . 

A color converting method for converting into 
a full color image at higher speed and lower cost after the 
image is compressed in a JPEG format is disclosed (for 
example, refer to IDS or Cross -Reference 3) . In this color 
converting method, the input color image is separated into 
a brightness element and a chrominance element, a color 
conversion number is obtained by using the chrominance 
element in which the brightness element is input and an 
output color space value output, a color conversion is 
conducted by applying the brightness element to a color 



converting f unction , and then an input color in a given color 
space is converted into an output color in another color 
space. 

The IDS or Cross-References 1 and 2 discloses 
the reversible color conversion by the integer operation , 
in the reversible color conversion, an error is not 
accumulated even if a conversion and a reversed conversion 
are repeated, and then the color conversion is reversibly 
conducted at a predictable accuracy. Therefore , it is just 
a simple data conversion. 

Moreover^ in a case of IDS or Cross -References 
3, the color converting process can be conducted at high 
speed, and an operation amount and cost can be reduced in 
the color converting process. Accordingly, a conversion 
accuracy can not be improved. 

That is, in JPEG and JPEG 2000, as describe 
above, an image subject to be compressed is processed by 
YCbCr (Y: brightness, Cb: Black color difference, Cr:Red 
color difference) data. In addition, in generally, the 
image display and a print developing process are processed 
by RGB and YMC (or YMCK) . Accordingly, compatibility for 
both RGB and YMC can be problem. In particular, in JPEG 
2000, when the electronic watermark for preventing form 
munipulaton and the electronic watermark for encryption are 
embodied, it is required to loss any data even if the 
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conversion and the reversible conversion are conducted / and 
maintain the compatibility of image data conversion. In 
particular, in JPEG 2000 where lessless compression and 
extension are available , RGB and YCbCr color conversion is 
a not-reversible conversion. It is required for a basic 
color conversion disclosed in IDS or Cross-References 1 and 
2 to guarantee the image quality. 

In addition, not only the conversion for the 
color space, a problem concerning a length conversion of 
the image data can be raised. For example, in general, data 
in a personal computer is process in a mm unit system but 
data in a printer is processed in a inch unit system. 
Accordingly , if the reversibility in the conversion and the 
reversed conversion is displaced, a page location may be 
displaced. 

[IDS or Cross-Reference 1] 

•Japanese Laid-Open Patent Application No. 9-6952 
[IDS or Cross -Reference 2] 

Japanese Laid-Open Patent Application No. 11-219428 
[IDS or Cross -Reference 3] 

Japanese Laid-open Patent Application No. 
2000-175061 

SUMMARY OF THE INVENTION 

It is a general object of the present invention 
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to provide a method for revozsibly transforming data format; 
an Image processing apparatus, a program for reversibly 
transforming data format , and a computer-readable 
recording medium thereof, in which the above-mentioned 
problems are eliminated. 

A more specific object of the present invention 
is to provide a a method for reversibly transforming data 
format, an image processing apparatus, a program for 
reversibly transforming data format, and a 
computer-readable recording medium thereof , in which a 
color space transforming process is conducted in a system 
where a lossless/lossy compression/extraction is conducted 
with respect to the image data in a color space. 

The above objects of the present invention are 
achieved by a method for reversibly converting a data format 
in that a forward transformation and a backward 
transformation are reciprocally conducted for data between 
unit systems having different resolution levels, wherein 
in the forward transformation and the backward 
transformation; a first unit system having a lower 
resolution level is used as a common unit system, and a 
reversible data conversion is conducted by an integer 
operation for data in the first unit system having the lower 
resolution level and data in a second unit system having 
a higher resolution level higher than the first unit system. 
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In the method, the first unit system is a first 
color space having the lower resolution level, and the 
second unit system is a second color space having the higher 
resolution level , wherein in a case in that a digital color 
conversion by quantizing the first color space having the 
lower resolution level and being analog, by using the first 
color space as the common unit system, the reversible data 
conversion is conducted by the integer operation, so that 
the backward transformation from first data in the first 
color space to second data in the second color space and 
the forward transformation from the second data in the 
second color space to third data in the first color space 
are conducted and the first data corresponds to the third 
data. 

In the method, the first data in the first color 
space are YCbCr data, and the second data in the second color 
space are RGB data, and the third data in the first color 
space are Y'Cb'Cr' , and in a case of converting to the 
Y' Co' Cr' data after the YCbCr data is converted into the 
RGB data, when the forward transformation is conducted from 
the RGB data to the Y'Cb'CR' data, a first formula is 
conducted by an integer operation of a color converting 
functions where the first formula is defined as 
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Y I 219x(299xR+587xG+114xB) + 16 x 255x1000 + 255x1000/2 I 
L 255x1000 J 

,■■ I 224 x 564 x (-299 x R - 587 x O + 886 x B) + 128 x 255 x 1000 x 1000 + 255 x 1000 x 1000/2 

CD = I — ■ — - . . 

L 255x1000x1000 
Cr = | 224x713x(701xR-587xG-114xB) + 128x25SxlO00xlO0O + 25SxlOOOxlOOO/2 
1 [ 255x1000x1000 

[ J : round a fractional part and the same in the following 

and when the backward transformation is conducted from the 
RGB data to the YCbCr data, a second formula is conducted 
by an integer operation of a color converting functions 
where the second formula is defined as 

„ I [219 x 1000 x (Cr -128)+ 713 x 224 x (Y -I6)]x 255 + 713 x 224 x 219/2 



G = 



713 x 224 x 219 
'713 x 224 x 587 x 564(Y - 16) 
- 299 x 219 x 564 x 1000 x (Cr - 128) 
-114 x 219 x 713 x 1000 x (Cb - 128) 



x 255 + 219 x 713 x 224 x 587 x 564/2 



219 x 713 x 224 x 587 x 564 



B I [219 x 1000 x (Cb - 128)-+ 564 x 224 x (Y - 16)]x 255 + 564 x 224 x 219/2 
I 564 x 224 x 219 



In the method, the first unit system is a first 
color space having the lower resolution level and the second 
unit system is a second color space having the higher 
resolution level, wherein in a case in that a digital color 
conversion by quantizing the first color space having the 
lower resolution level and being analog, by using the first 
color space as the common unit system, the reversible data 
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conversion is conducted by the integer operation, so that 
the backward transformation from first data in the first 
color space to second data in the second color space and 
the forward trans formation from the second data in the 
second color apace to third data in the first color space 
are conducted and the first data corresponds to the third 
data. 

In the method, the first data in the first color 
space are YCbCr data , and the second data in the second color 
space are RGB data, and the third data in the first color 
space are Y'Cb'Cr', and in a case of converting to the 
Y'Cb'Cr' data after the YCbCr data is converted into the 
RGB data, the forward transformation from the second data 
in the second color space to the third data in the first 
color space is conducted by an integer operation where MAX 
denotes a maximum value of the resolution level and Y, M, 
and C are set to be Y=MAX-B , M=MAX-G, and C=MAX-R in the 
first formula, and the backward transformation from the 
first data in the first color space to the second data in 
the second color space is conducted by an integer operation 
where MAX denotes a maximum value of the resolution level 
and Y, M, and C are set to be Y=MAX-B, M=MAX-G , and C=MAX-R 
in the second formula. 

In the method, the first unit system is a first 
color space having the lower resolution level and the second 
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unit system is a second color* space having the higher 
resolution level, wherein in a case of conducing a color 
conversion in accordance with an international standard in 
which the data format for converting an analog video signal 
into digital data is specified, by using the first color 
space as the common unit system, the reversible data 
conversion is conducted by the integer operation, so that 
the backward transformation from first data in the first 
color space to second data in the second color space and 
the forward transformation from the second data in the 
second color space to third data in the first color space 
are conducted and the first data corresponds to the third 
data. 

In the method, the first data in the first color 
space are YCbCr data, and the second data in the second color 
space are quantised R(d)G(d)B<d) data, and the third data 
in the first color space are Y'Cb'Cr', and in a case of 
converting to the Y'Cb'Cr' data after the YCbCr data is 
converted into the RGB data, in a conversion from an analog 
R(a)G(a)B(a) data to the YCbCr data, the forward 
transformation from the YCbCr data by conducting a color 
converting function in a third formula defined as 
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J 219 x (77 x R(a) + 150 x G(a) + 29 x B(a)) + 16 x 256 x 256 + 256 x 128 I 
I 256 x 256 J 

_ I 2 19 x (- 44 x R(a) - 87 x G(a) + 1 3 1 x B(a)) + 1 28 x 256 x 256 +. 256 x 128 I 
I 256 x 256 J 

I 219x(l31xR(a)-U0xO(a)-21xB(a)) + 128x256x256 + 256xl28 I 
r "[ 256x256 J 

to the R(d)G(d)B (d) data is conducted by an integer 
operation of a color converting function in a fourth formula 
defined as 



R(d)-[i 
Q(d)-[J 
B(d)=|l 



(16772821 x Y + 22904709 x Cr - 41320 x Cb - 2926513792) x 2 + 16772821 

16772821x2 

,U _| (470873 x Y - 329527 x Cr - 157064 x Cb + 62283648)x 2 + 470873 
* ' l 470873x2 

| (16772821 x Y - 102267 xCr + 29047960 x Cb - 3705048704)x 2 + 16772821 

16772821x2 



and when the backward transformation is conducted from the 
R<d)6(d)B(d) data to the Y'Cb'Cr' data, a fifth formula is 
conducted by an integer operation of a color converting 
functions where the second formula is defined as 



Y' = 



77 xR(d)+150 xG(d)+29 xB(d)+128 
256 

■ I -44x R(d)-87 xG(d)+l31 xB(d)+128 x 256 +128 I 
= L 256 J 

, I 131 x R(d)-110 x G(d)- 21 x B(d) + 128 x 256 + 128 I 
r I 256 J . 

In the method, the first unit system is a first 
color space having the lower resolution level, and the 
second unit system is a second color space having the higher 
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resolution level, wherein in a case of conducting an 
original color based on the brightness and the color 
difference, by using the first color space as the common 
unit system/ the reversible data conversion is conducted 
by the integer operation, so that the backward 
transformation from first data in the first color space to 
second data in the second color space and the forward 
transformation from the second data in the second color 
space to third data in the first color space are conducted 
and the first data corresponds to the third data. 

In the method, the first data in the first color 
space are Y (o) Cb (o) Cr (o) data, and the second data in the 
second color space are R(o) G (o)B (o) data, and the third data 
in the first color space are Y (o) ' Cb (o) ' Cr (o) ' , . and in a 
case of converting to the Y (o) ' Cb (o) ' Cr (o) ' data after the 
Y(o)Cb(o) Cr<o) data is converted into the R (o) G (o) B (o) data, 
when the forward transformation is conducted from the RGB 
data to the Y (o) ' Cb (o) ' CR (o) ' data, a sixth formula is 
conducted by an integer operation of a color converting 
functions where the sixth formula is defined as 
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y(o ) |^ 2x ^»' xR < ;o )- | -( p - x »- y ^) xG ( o ) + y» xB ^))- f - p j 

j x 2x(D — y M )— x M x R(o)- (D - x M - y M )x G(o) + (b - y M > (B (o) + 1) 



Cb(o) = 



Cr(o) = 



2x(D- y ;j 
| MAX „, B +l j 

| M^BS&llj x 2x (D - x M )+ (D -x»)x (r(o)+1)-(D - x„ - y „ )x G(o)- y „ x B(o) 



2*(D-*m) 



and when the backward transformation is conducted from the 
R(o)G(o)B(o) data to the Y(o)Cb<o)Cr (o) data, a seventh 
formula is conducted by an integer operation of a color 
converting functions where the seventh formula is defined 



K(Q-[ 2x ^ xY(0)+2X 2 (P D " XM)xCr(0))+D J 

r2x((D-x M -y M )xDxY(o)-2xy M x(D-y M )xCb(o)-2xx M x(D- XM )xCr(o))^ 
{ -t-(D-x M -y M )xD ) 



G(o) 



2x(D-x M -y M )xD 



2x(PxY(o)+2x(P-y M )xCb(o))-i-P 
2xP 



J 



In the method as claimed in claim 9, wherein 
a color conversion for the original color based on the 
brightness and the color difference, the first data in the 
first color space are Y (o) Cb (o) Cr (o) data, and the second 
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data in the second color space are R(o)G(o)B(o) data, and 
the third data in the first color space are 
Y(o) 'Cb(o) 'Cr(o) ' , and in a case of converting to the 
Y(o) 'Cb(o) 'Cr(o) ' data after the Y (o) Cb (o) Cr (o) data is 
converted into the R<o)G(o)B(o) data, when the forward 
transformation is conducted from the RGB data to the 
Y (o) ' Cb(o) ' CR(o) ' data, an eighth formula is conducted by 
an integer operation of a color converting functions where 
the eighth formula is defined as 

, v I (299 xR(o)+587 xG(o)+U4 xB(o))+500 
^'l 1000 
Cb (o) - 1^ 128 x2x886 -299 xR(o)-587 x G(o)+ 886 x (B(o)+ l) j _ 12g 

( v I 128 x 2 x 701 + 701 x (R(o)+ l)- 587 x G(o)- 114 x B(o) I - 
t[9) ~[ 2x701 J 

and when the backward transformation is conducted from the 
R(o)G(o)B(o) data to the Y (o) Cb (o) Cr (o) data, a ninth 
formula is conducted by an integer operation of a color 
converting functions where the ninth formula is defined as 

_/* I (l 000 x Y(o) + 1402 x Cr(o)) + 500 

R(oj l looo . 

q/ \ J (587 x 1000 x Y(o) - 2xll4x886x Cb(o) - 2 x 299 x 701 x Cr(o)) + 587 x 500 
*°'~|_ 587x1000 J 

n I (lOOOx Y(o)+1772xCb(o))+500 I 
*°'"|. 1000 J 



In the method, maximum values and minimum value 
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are limited in available ranges for the R(o)G(o)B(o) data 
and the Y (o) Cb (o) Cr (o) data. 

In the method, wherein the first unit system 
is a BMU unit system using an inch unit system and the second 
unit system is a 1/100 mm unit system using a meter unit 
system, and by using the first color space as the common 
unit system, the reversible data conversion is conducted 
by the integer operation, so that the backward 
transformation from first data in the first color space to 
second data in the second color space and the forward 
transformation from the second data in the second color 
space to third data in the first color space are conducted 
and the first data corresponds to the third data. 

In the method, after the first data shown in 
the BNU unit system is converted into the second data shown 
in the 1/100 mm using the meter unit system, a tenth formula 
is conducted to convert to the third data showing the BMU 
unit system where the tenth formula is defined as 
<l/100mm data> =[(2540x<BMU data> + 600)/1200 J 

<BMUdata>=L(l200x<l/100mradata> + 1270)/2540 J 

In the method, the integer operation conducts 
the reversible conversion using powers of 2. 

In the method, the forward transformation from 
the R(o)6(o)B(o) data to the Y (o) ' Cb (o) ' Cr (o) ' data is 
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conducted by an eleventh formula using powers of 2 derived 
from the sixth formula where the eleventh formula is defined 
as 



(l28 x 2 U 


2 12 

)- 691 xR(o)- 1357 xG(o)+2" 


xB(o)+2" 


(128 x 2 12 ) 


|-2 n x R(o)-1715 xG(o)-333 


xB(o)+2" . 




.2" 





Cbi 
Cr( 

and the backward transformation from the Y (o) Cb (o) Cr (6) 
data to the R(o)G(o)B(o) data is conducted by a twelfth 
formula using powers of 2 derived from the seventh formula 
where the twelfth formula is defined as 

(l28 x2")+2 12 x Y(o2^5743xCr(o)+2 n m 
2 

■', . (l28x2 ,2 )+2 12 x Y(o)- 1410 xCb(o)- 2925 xCr(o)+2 n ' 
G(o)=^ 1- ^ -s- 128 

B(o)= ( 128 x2l2 ) +2 " x Y(o)+7258 xCb(o)+2" 12g 

In the method, the powers of 2 is conducted by 

bit shifts. 

In the method, the forward transformation from 
the R(o)6<o)B(o) data to the Y (o) ' Cb (o) ' Cr (o) ' data is 
conducted by a thirteenth formula using powers of 2 derived 
from the sixth formula where the thirteenth formula is 
defined as 
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Y(o)=(((128«12)+1225X R(o)+2404x G(o)+467 X B(<»)+(1«11))»12>-128; 
Cb(a)=(((128«12)-691 x R(©)-1357 x G(o)+204B x B(o)+(1«1 1))»12)-128: 
CKo)=(((12B«1 2)-2048 x R(«,)- 1 71 5 x G(e>)-333 x B(e)+(1 «1 f ))»12)-1 2B; 

and the backward transformation from the Y (o) Cb (o) Cr (o) 
data to the R (o) G (o) B (©) data is conducted by a fourteenth 
formula using powers of 2 derived from the seventh formula 
where the fourteenth formula is defined as 

R(o)= (((l28«12) + 4096 x Y(o)+ 5743 x Cr(o)+ 1)>>12)- 128; 

G(o)= (((l28«12)+4096 x Y(o)- 1410 x Cb(o)-2925 x Cr(o) + (l«ll))»12)-128; 

B(o)= (((l28<<12) + 4096 x Y(o)+ 7258 x Cb(o)+ (l«l l))»12)- 128; 

« 

According -to the present invention, in the 
forward transformation and the backward transformation , a 
unit system having a lower resolution revel is used as the 
common unit system and the data conversion is conducted by 
the integer operation between unit systems having different 
resolution levels , by using the common unit system. 
Therefore, a unit system having a higher resolution level 
is sufficiently processed and then it is possible to realize 
the reversible conversion where data can be perfectly 
restored. 

In particular , regarding the digital color 
conversion using a quantization of analog Y, Cb, Cr signals 
between the unit system for the color space RGB and the unit 
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system for the color space YCbCr, the unit system for the 
color space YCbCr having a lower resolution level is used 
as the common unit system. Thus, since the unit system for 
the color space RGB is sufficiently processed, it is 
possible to realize the reversible conversion in a color 
space where data can be perfectly restored. 

Moreover, it is possible to realize the 
reversible conversion in the color space where data can be 
perfectly restored by a predetermined color converting 
function . 

Furthermore, in a case in that a unit system 
for a print in such as a color space YMC is used, it is 
possible to realize the reversible conversion in the color 
space where data can be perfectly restored. 

In particular, regarding the color conversion 
in accordance with the international standard specifying 
a data format for converting an analog video signal into 
the digital data between the unit system in the color space 
RGB and the unit system in the color space YCbCr, the color 
space YCbCr having a lower resolution level is used as the 
common unit system. Therefore, since the unit system in 
the color space RGB having a higher resolution level can 
be sufficiently processed, it is possible to realize the 
reversible conversion of the color space where data can be 
perfectly restored. 
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In particular, regarding the color conversion 
of the original color based on the brightness and the color 
difference between the unit system in the color space RGB 
and the unit system in the color space YCbCr , the color space 
YCbCr having a lower resolution level is used as the common 
unit system. Therefore , since the color space RGB having 
a higher resolution level can be sufficiently processed, 
it is possible to realize the reversible conversion in the 
color space where data can be perfectly restored. 

Moreover, in the available ranges for 
R(o)6(o)B(o) data and ?(o) Gb(o)Cr(o) data, a limited range 
is provided in that a maximum value and a minimum value are 
limited. Therefore, without being out of the range of the 
cplor converting function, it is possible to conduct the 
compatible conversion where data can be perfectly restored. 

In particular, regarding the length conversion 
between the BMU type data using an inch unit system and the 
1/100 mm type data using a meter unit system in the personal 
computer and the printer, the inch unit system for the BMU 
type data having a lower resolution level is used as the 
common unit system. Therefore , since the meter unit system 
for the 1/100 mm type data having a higher resolution level 
can be sufficiently processed, it is possible to realize 
the reversible conversion where length data can be 
perfectly restored. 
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The above objects of the present: invention can 
be achieved by an image processing apparatus that realizes 
-the process in the method for reversibly transforming data 
format, by a program code for causing a computer to conduct 
processes described above in the method or by a 
computer-readable recording medium recorded with the 
program code. The image processing apparatus is not 
limited to a single function device such as a digital camera, 
a digital copier , and an MFP (Multi -Functional Printer) . 
Alternatively, the image processing apparatus can have a 
system configuration of a server /client system or a like . 

BRIEF DESCRIPTION OF THE DRAWINGS 

In the following; embodiments of the present 
invention will be described with reference to the 
accompanying drawings. 

FIG.l is a diagram schematically showing a 
system configuration according to an embodiment of the 
present invention; 

FIG, 2 is a diagram schematically illustrating 
an application to JPEG 2000; 

FIG.3 is a diagram schematically illustrating 
an application to a usage requiring a reversible color 
conversion; 

FIG. 4 is a diagram schematically illustrating 
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results of an analog type color conversion; 

FIG * 5 is a diagram schematically illustrating 
results of a digital type color conversion; and 

FIG. 6 is a diagram schematically showing 
results of the color conversion applied to JPEG 2000. 

DESCRIPTION OF THE PREFERRED EMBODIMENT 

In the following, an embodiment of the present 
invention according to a document processing apparatus will 
be described with reference to the accompanying drawings . 
The present invention can be applied to various image 
processing systems and devices such as a color television 
broadcast, a digital video, a digital camera, a digital 
copier, an MFP (Multi -Functional Printer) f various image 
processing application software for a personal computer, 
a printer, a scanner, and a like. For the sake of 
convenience, just the most essential part of the present 
invention is explained and shown in the accompanying 
drawings regarding the various image processing systems and 
image processing devices. In addition/ a data format 
reversible converting method according to the present 
invention can be realised by a program which can be 
executable in these various image processing system and 
image processing devices, and various computer -readable 
recording media recording the program. 
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[System Configuration] 

The present invention can be included in a 
compressing/ extending system having a coding part and a 
decoding part that are provided to the image processing 
apparatus (any one of the various image processing system 
and image processing devices) . PIG.l illustrates a system 
configuration according to the present invention . A coding 
part 1 encodes input data to generate compressed data 
(encoded data) , and a decoding part 2 decodes the encoded 
data to generate original input data. In the coding part 
1, a color space transformation 4 is conducted to convert 
a first color coordinate system into a second color 
coordinate system prior to a compression 5. Similarly, in 
the decoding part 2, a color space transformation 7 is 
conducted from the second coordinate system to the first 
coordinate system after an extension 6. 

The input data for the compressing/expressing 
system may be an image (static image or dynamic image) , 
graphics data, and data in other various data formats. In 
this case, the input data is digital signal data. 
Alternatively, data in other formats can be the input data. 
A source of the data may be based on a communication path 
or a memory for the coding part 1 and/or the decoding part 
2. This system can be configured as a lossless 
compressing/ extending system or to conduct a lossy 



26 



compression/extension. 

In addition, aa described above, configuration 
factors of the coding part 1 and/or the decoding part 2 may 
be realized by hardware, software, or a combination of both 
hardware and software. The color space transformations 4 
and 7 are realized as a data format lossless transforming 
part or that function, as follows. 

[Application for Analog Type Compatible Color 
Transformation] 

An example of an application to a digital color 
transformation by quantizing analog Y, Cb, and Cr signals 
(Y : brightness , Cb:B (Black) color difference, and Cr:R 
(Red) color difference) will be described according to an 
embodiment of the present invention. In this case, 
regarding the digital color transformation by quantizing 
the analog Y, Cb, and Cr signals between an unit system in 
a color space RGB and an unit system in a color space YCbCr, 
the unit system of the color space YCbCr having a lower 
resolution is used as a common unit system to maintain a 
compatibility of an integer operation with the color space 
RGB. That is, between the unit system of the color space 
YCbCr having a lower resolution and the unit system of the 
color space RGB used for a display by using the unit system 
of the color space YCbCr, the following data conversion is 
conducted in an order of a backward transformation and a 
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forward transformation: 

YCbCr data -> RGB data -> Y'Cb'Cr data. 
The data conversion is reversibly conducted by the Integer 
operation so as to satisfy Y=Y' , Cb=Cb' , and Cr=Cr' . The 
inventors found a color converting function for conducting 
the forward transformation and another color converting 
function for conducting the backward transformation in the 
data conversion for conducting the backward transformation 
and the forward transformation to realize a YCbCr data -> 
RGB data -> Y'Cb'Cr data conversion. It should be noted 
that each of R, G , and B element values used in the embodiment 
is independent from each other and is one of values from 
2 to 253 (expressed in 8 bits) . 

First, as an RGB data -> YCbCr data conversion, 
the following color converting function can be applied to 
a digital color converting function quantizing the analog 
Y, Cb, and Cr signals cited in ITU-R BT.601 being an 
international standard that specifies a data format when 
an analog video signal specified for the color television 
broadcast or the like is converted into a digital data. 
Regarding the digital color converting function , a 
converting process is conventionally conducted by a real 
number operation shown in the following formula 15: 
[Formula 15] 



28 



Y=219x[o.299x — + 0.587x^ + 0.114*^ | + 16 
[ 255 255 255 J 

Cb = 224 x|a564xf- 0.299 x — -0.587 x^ + 0.886 x^|| + 128 
[ ^ 255 255 255 J J 

Cr = 224 x |^0.713 x (a701 x ^ - 0.587 x ^ - 0.1 14 x i j j + 128 

However, when the formula 15 is calculated by a computer, 
errors occur due to an overflow o£ an arithmetic operation. 
On the contrary, in this embodiment, for the RGB data -> 
VCbCr data conversion, instead of using the formula 15, the 
following formula 16 is applied so as to conduct the RGB 
data -> XCbCr data conversion by an integer operation: 
[Formula 16] 

| 219x(299xR^587xG+114xB)+16 x 255xl000 ^255xlOOO/2 I 
Y "[ ~ 255x1000 " j 

I 224x564x(-299xR-587xG+886xB) + 128x255xl000xl0(K) + 255x1000x1000/2 I 
L 255x1000x1000 J 

I 224x713x(70lxR-587xG-114xB)+128x255xlO00xlOO0 + 255xl0O0xl000/2 I 
~[ 255x1000x1000 J 

[ J : roun<J a fractional part and the same in the following. 

In the formula 16, an addition operation and 
a multiplication operation, which do not cause an error 
unless the overflow occurs , are conducted in an earlier step . 
A division operation, which easily cause the overflow and 
the error , is conducted in a later s tep . In addition , after 
half denominator is added, a fractional part is rounded, 
so as to count fractions over 1/2 as one and disregarding 
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the rest. 

That is, in a Y data converting expression in 
the formula 16, after "255x1000/2" showing half the 
denominator is added to a numerator, the numerator is 
divided by the denominator "255x1000" and then a fractional 
part is rounded. In Cb and Cr data converting expressions 
in the formula 16, "255x1000/2" showing half a denominator 
is added to a numerator, a result from calculating the 
numerator is divided by the denominator "255x1000x1000", 
and then a fractional part is rounded. 

It is possible to conduct an arithmetic 
operation in accordance with the formula 16 by integers and 
to obtain an operation result by an integer value. 

On the other hand, the digital color converting 
function, which is cited in the international standard 
XTU-R BT . 601 and quantizes the analog X, Ob, and Cr signals , 
can be applied to a YCbCr data <- RGB data conversion. 
Conventionally, a converting process by a real number 
operation shown in the following formula 17 is used: 
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[Formula 17] 

Cr-128 Y-16' 



;255 



713 x 224 219 _ 
Y-16 0.299 x(Cr- 128) 0.114x(Cb-128) 1 2g5 
219 0.713 x 224 x 0.587 ~ 0.564 x 224 x 0.587 J* 



Cb-IM 
'0.564 x 224 219 



However, when the formula 17 is calculated by a computer , 
errors occur due to an overflow of an arithmetic operation. 
On the contrary, in this embodiment, for the YCbCr data -> 
RGB data conversion, instead of using the formula 17, the 
following formula 18 is applied so as to conduct the YCbCr 
data ~> RGB data conversion by an integer operation: 
[Formula 18] 



4 



G = 



[219 x 1000 x (Cr - 128)+ 713 x 224 x (Y - 16)]x 255 + 7 13 x 224 x 219/2 

713 x 224 x 219 

7 13 x 224 x 587 x 564(Y - 16) 
- 299 x 2 19 x 564 x 1000 x (Cr - 128) 
-114x219x713xlOO0x(Cb-128) 



x 255 + 219 x 713 x 224 x 587 x 564/2 



219 x 713 x 224 x 587 x 564 



B I [219xl000x(Cb~128) + 564 x 224x(Y-16)]x255 + 564 x 224 x 219/2 I 
~|_ 564 x 224 x 219 J 

In the formula 18, an addition operation and 
a multiplication operation/ which do not cause an error 
unless the overflow occurs, are conducted in an earlier step . 
A division operation, which easily cause the overflow and 
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the error, is conducted in a later step. In addition, after 
half denominator is added, a fractional part is rounded, 
so as to count fractions over 1/2 as one and disregarding 
the rest. 

That is, in a R data converting expression in 
the formula 18, after "713x224x219/2" showing half the 
denominator is added to a numerator, the numerator is 
divided by the denominator "713x224x219" and then a 
fractional part is rounded. In a G data converting 
expression in the formula 18, "219x713x224x587x564/2" 
showing half a denominator is added to a numerator, a result 
from calculating the numerator is divided by the 
denominator "219x713x224x587x564" , and then a fractional 
part is rounded. In a B data converting expression in the 
formula 18, "564x224x219/2" showing half a denominator is 
added to a numerator, a result from calculating the 
numerator is divided by the denominator "564x224x219", and 
then a fractional part is rounded. 

It is possible to conduct an arithmetic 
operation in accordance with the formula 18 by integers and 
to obtain an operation result by an integer value . 

As described above, in the YCbCr data -> RGB 
data -> Y'Cb'Cr' data conversion, by using the color space 
YCbCr as the common unit system, the backward 
transformation at the YCbCr data -> RGB data side is 
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conducted by using the color converting expression shown 
by the formula 16, and the forward transformation at the 
RGB data -> Y'Cb'Cr' data side is conducted by using the 
color converting expression shown by the formula 18 . Table 
1 shows a color conversion results of the backward 
transformation by the formula 16 and the forward 
-transformation by the formula 18. 
[Table 1] 



example 1. 



YCbCr 
data 



16 (Y) 
128 (Cb) 
128 (Cr) 



(by formula 1 8) 



RGB 
data 

0(G) 
0(B) 
0(R) 



(by formula 16) 



restored 

YCbCr 

data 

16 (Y) 
128 (Cb) 
128 (Cr) 



example 2. 



126 (Y) 

12B(Cb) 

128(Cr) 



128(G) 
128(B) 
128 (R) 



126 (Y) 

128(Cb) 

l28(Cr) 



example 3. 



107(Y) 
191 (Cb) 
208 (Cr) 



16(G) 

233 (B) 

234 (R) 



107(Y) 
191 (Cb) 
208 (Cr) 



example 4. 



example 5. 



112(Y) 
139 (Cb) 
123(Cr) 
218 (Y) 
128 (Cb) 
128 (Cr) 



112(G) 
134(B) 
104(R) 
235 (G) 
235 (B) 
235 (R) 



112(Y) 
139 (Cb) 
123 (Cr) 

218 (Y) 
128(Cb) 
128 (Cr) 



According to the color conversion result shown 
in Table 1, Y=Y' , Cb=Cb' , and Cr=Cr' are maintained. Thus , 
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it can be seen that a data compatibility (reversibility) 
concerning a color conversion is guaranteed, 

lExample of Application to Color Space YMC Used for 

Print] 

In the previous example , the application to the 
unit system of the color space RGB used for a display is 
described. Also, the present invention can be applied to 
a unit system of the color space YMC used for a print. In 
this case, Y, M, and C are set to be Y^MAX-B , M=*dAX-G, and 
C=max-R in the formula 40 and the formula 42 , so as to conduct 
the integer operation. MAX denotes a maximum value of the 
resolution level . For example , it the maximum level is 256 , 
MAX is set to be MAX=256 . 

[Example, of Application to Digital Type 
Compatible Color Conversion] 

Another application to a color conversion 
formula cited in the international standard ITO-R BT . 601 
as an example. In this case, the unit system of the color 
space YCbCr is used as the common unit system. The data 
conversion of the backward transformation and the forward 
transformation shown as 

YCbCr data -> R(d)G(d)B(d) data -> Y'Cb'Cr' data 
is conducted between the unit system of the color space YCbCr 
and the unit system of the color space RGB used for the 
display, to conduct the reversible data conversion by the 
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integer operation so that Y=Y' , Cb=Cb' , and Cr=Cr' are 
satisfied/ In the data conversion such as the YCbCr data 
-> R(d)G(d)B<d) data -> Y'Cb'Cr' data conversion, the 
inventors found a color converting function for conducting 
an R(d) 6(d) 8(d) data ~> YCbCr data conversion, a color 
converting function for conducting a forward 
transformation of a YCbCr data <- R(d)6(d)B(d) data 
conversion, and a color converting function for conducting 
a backward transf omation of an R(d)G(d)B(d) data -> 
Y'Cb'Cr' data conversion. 

An original color converting function is a 
color converting function cited in the international 
standard ITU-R BT.601, In the following, conventional 
converting formulae will be described. 

a. converting formula from analog RGB (R(a) , G(a) , 
B(a)) signals to quantized digital (R(d) , G(d) , B(d)) 
signals 

where 0^ R(a) S 256, 0^ G(a) ^ 256 , 0^ B(a) ^ 

256, 

[formula 19] 

R(d)=219 * R(a)/256 + 16 
G(d)=219*G(a)/256 + 16 
B(d)=21 9 * B(a)/256 + 1 6 

b. converting formula from the quantized digital 
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(R(d) , 6(d), B(d)) signals to digital Y, Cb, Cr signals 
[formula 20] 

Y=(77 * R(d)/256)+(1 50 * G(d)/256)+(29 * B(d)/256) 
Cb=-(44 * R(d)/256)-(B7 * G(d)/256)+(131 * B(d)/256)+128 
Cr=(131 * R(d)/256)-(110* G(d)/256)-(21 * B(d)/256)+128 

Conventionally, the above formula 20 is used for a 
converting process. On the other hand, according to the 
present invention, the R(a)G(a)B(a) data -> YCbCr data 
conversion is conducted by an operation in accordance with 
a color converting function shown by the following formula 
21: 

[Formula 21] 

y I 219x(77xR(a)-t-150xG(a)+29xB(a)) + I6x2S6x2S6 + 256xl28 | 
L ~" " 256x 256 J 

rb _ I 21 9 * (~ 44 * Rfr) ~ 87 * G(a)+ 131 x B(a))+ 1 28 x 256 x 256 + 256 x 128 I 
L ™ 256 x 256 J 

Cr = | 219x (131 xR(a)-U0xG(a)-21xB(a))-f 128 x 256 x 256 + 256x128 I 
~L 256x256 J 

The YCbCr data -> R(d) G (d) B <d) data conversion is conducted 
by an integer operation in accordance with a color 
converting function shown by the following formula 22 : 
[Formula 22] 
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x | (l6772821xY + 22904709xCr-41320xCb-29265l3792)x2 + 16772821 
R W"[ 16772821x2 

I (470873 x Y - 3295 27 x Cr - 157064 xCb + 62283648)x 2 + 470873 I 
G W = [ ; 470873 x 2 ~~J 

v I (l6772821xY-102267xCr + 29047960xCb-3705048704)x2 + 16772821 I 
B(d) = | — ■ 1677 2821x2 J 



The R(d)G(d)B(d) data -> Y'Cb'Cr' data conversion is 
conducted by an integer operation in accordance with a color 
converting function shown by the following formula 23: 
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[Formula 23] 

Y 



I 77 x R(d) + 1 50 x G(d) + 29 x B(d) + 128 I 
cb , |^-44x R(d) - 87 x G(d) 4- 1 3 1 x B(d) + 128 x 256 + 1 28 j 
Cft | l31x R(d)-110x G(d)~ 21 x B(d)+128x 256 + 128 j 



Accordingly, by using the color space YCbCr as 
the common unit system, for the data conversion of the 
backward transformation and the forward transformation of 
the YCbCr data -> R(d)6(d)B(d) data -> Y'Cb'Cr' data 
conversion, the R(a)G(a)B(a) data -> YCbCr data conversion 
is conducted by using the color converting function shown 
by the formula 21, the YCbCr data -> R(d)G(d)B(d) data 
conversion is conducted by using the color converting 
function shown by the formula 22, and the R(d)G(d)B(d) data 
-> Y'Cb'Cr' data conversion is conducted by using the color 
converting function shown by the formula 23 . Table 2 shows 
color conversion results from the formulae 21, 22, and 23. 
[Table 2] 
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(by formula 22) (by formula 23) 



example 1. 

example 2. 

example 3. 

example 4. 
example 5. 

in the Table 2 , the YCbCr data are restored. That is , Y=Y' , 
Cb=Cb' , and Cr=Cr' are maintained so that the data 
compatibility (reversibility) concerning the color 
conversion can be guaranteed. 

[Application to Original color converting function 
Bing Based on Brightness and Color Difference : JPEG 2000] 

In a JPEG 2000 algorithm, a color space 
converting process is conducted by converting the RGB data 
into the YCbCr data being based on the brightness and color 



YCbCr > R(d)G(d)B(d) _ ^ 

data data datg 

16 (Y) 16(G) 16 (Y) 

128(Cb) 16(B) 128 (Cb) 

128(Cr) 16 (R) 128 (Cr) 



128 00 128(G) 128 (Y) 

128(Cb) 128(B) 128 (Cb) 

128(Cr) 128(R) 128 (Cr) 



106(Y) 16(G) 106 (Y) 

202 (Cb) 234(B) 202 (Cb) 

222 (Cr) 234 (R) 222 (Cr) 



112(Y) 112(G) 112(Y) 

139 (Cb) 131(B) 139 (Cb) 

123 (Cr) 105(R) 123 (Cr) 

235 (Y) 235(G) 235 (Y) 

128(Cb) 235(B) 128(Cb) 

128 (Cr) 235 (R) 128(Cr) 

According to the color conversion results shown 
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difference or conducting the backward transformation . 
However, this color space converting process can not be 
reversible as described in Description of the Related Art. 

A. That is, as the color converting function of the 
conventional JPEG algorithm, 

a. conversion from RGB (R(0) , 6(0) , B(0) ) signals to 
YCbCr (Y(o), Cb(o) , Cr(o)) signals 

The following converting formula shown in 
formula 24 is used for the forwarding transformation: 
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[Formula 24] 









Cb(o) 




G(o) 


Cr(o)_ 




B(o). 



A matrix denoted by a letter A in the formula 2 4 Is described 
as follows : 



[Formula 25] 



W-y) 

2 



1-x-y 


y 


1-x-y 


i 


2(1 -y) 


2 


1-x-y 


y 


2(l-x) 


2(l-x)J 



.A J = 



2(l-y) 



2(1- x) 



! 2y(l-y) 2x(l-x) 
1-x-y 1-x-y 



|A| = . 



1-x-y 



*(i-*)(i-y) 

In the formula 25, x denotes a weighted coefficient of a 
red color, y denotes a weighted coefficient of a blue color, 
and 1-x-y denotes a weighted coefficient of a green color 
(where 0 Sx<l . 0 ^y<l , and x+y<l ) . 

b. conversion from YCbCr (7(0) , Cb(0) , Cr(0)) signals 
to RGB (R(0), 6(0), B(0)) signals 

The following formula 26 is used for the 
backward transformation: 
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[Formula 26] 







[y(q)] 


G(o) 


=A J 


Cb(o) 


B(o). 




_Cr(o)_ 



A matrix denoted, fay a letter A in the formula 26 is the same 
matrix shown by the formula 25. 

The matrix A shown by the formula 25 will be 
described in detail. When x=299/1000=weighted 
coefficient of the red color, Y=114/1000=weighted 
coefficient of the blue color, and l-x-y=587/1000=weighted 
coefficient of the green color, the matrix A is expressed 
as follows: 
[Formula 27] 

1 0 1.402 



' 0.299 0.587 0.114] 

-0.16874 ^0.33126 0.5 
0.5 -0.41869 -0.08131 



1 -0.34414 -0.71414 
1 1.772 0 



The matrix A is a converting matrix conventionally used as 
non-reversible color conversion in the JPEG 2000. 

B. on the other hand, in the embodiment according to 
the present invention, similar to the previous applications , 
it is possible to conduct a reversible conversion by integer 
values for the color conversion in the JPEG 2000 algorithm. 
Thus , it is possible to obtain all variables and all 
coefficients and to obtain an operation result by an 
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integer. 

When x indicating the weighted coefficient of 
the red color, y indicating the weighted coefficient of the 
blue color, and 1-x-y indicating the weighted coefficient 
of the green color are expressed by x=^Xm/D and y=y M /D (where 
Xm, ym, and D are all integers, and an index M is simply 
a part of a variable name, for example, the index M is used 
to distinguish "x" from 'V) , x^, y M and D are expressed 
under limitations of x and y concerning the formula 25 
described above as follows : 
[Formula 28] 

O^XtfCD, O^Ym<D, and ^4+Ym<D. 

In this case , by using coefficients concerning 
a matrix, the color converting process is conducted as 
follows . It should be noted that " (o) " of R(o) , 6 (o) , B (o) , 
Y(o), Cb(o) , and Cr(o) is additionally provided to 
distinguish from the variable names for general R, G, B, 
Y, Cb, and Cr used in other converting formulae described 
above, in order for conversion data in this case according 
to the present invention to indicate "original color being 
based on the brightness and the color difference". Thus, 
"(°)" is simply a part of a variable name. 

a. conversion from RGB (R(0) , G(0) , B(0)) signals to 
YCbCr (Y(o), Cb(o), Cr(o>) signals 

A color converting function show in the 
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following formula 29 is conducted by the integer operation 
for the forward transformation. 
[Formula 29] 

y ^ ^ 2x(x M xR(o)+(P-x M -y M )xC(o)+y M xB(o))+D j 



Cb(o) = 



|^^^^j x 2x(D-y M )-x M xR(o)-(D-x M -y u )xG(o) + (D-y M )x(B(o) + l) 



Cr(o) = 



2x(D-y M ) 
2x(D-K M ) : _ _ 



When MAXrgb denotes a level number (integer) and 
a first limitation concerning input R(o.) , G(o) , B(o) data, 
which are used to realize reversibility and compatibility 
described later in the converting process is eliminated, 
R(o) , G(o), and B(o) are expressed as follows: 
[Formula 30] 

0^R(o) 0^G(o) ^MAXrcjb - 1 , 

O^B (o) gMAX RC B-X . 

For example, when x=29&/1000=weighted 
coefficient of the red color, y=114/1000=weighted 
coefficient of the blue color, l-x-y=587/l000=weighted 
coefficient of the green color, and levels in 

order to meet the JPEG algorithm, the formula 29 concerning 
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the RGB (R(o) , G (o) , B (o) ) signals -> YCbCr (Y(o) , Cb(o) , 
Cr(0)) signals conversion can be replaced with the 
following color converting function shown in formula 31 : 
[Formula 31] 

(299 x R(o) + 587 x G(o) + 1 14 x B(o)) + 500 | 



Y(o). 
Cb(o) = 
Cr(o) = 



1000 J 
128 x 2 x 886 - 299 x R(o) - 587 x G(o) + 886 x (B(o) + 1) I 

2x886 J" 
1 28 x 2 x 701 + 701 x (r(o) + 1) - 587 x G(o)- 1 14 x B(o) I 

2x701 



By combining the formula 31 and formula 38 
described later , it is possible to improve the color 
converting function, which is known as not-reversible in 
the JPEG 2000 algorithm to be a reversible color converting 
function. 

In addition, by applying a bit shift to the 
color converting function shown by the formula 31, it is 
possible to improve an operation speed. For example, since 
the level number of the RGB is 256 and an available range 
is from 0 to 255, the formula 31 can be expressed in C 
language where MAX_KAICH0^155 as a maximum level value , 
bitSHIFT=12 (the number of bit shifts) as exponent© of 
powers of 2 , iRORS, iGORS, and iBORS denote R(o> , G(o) , and 
B(o) , and iYORS, iCbRS, and iCrRS denote Y(o) , Cb(o) , and 
Cr(o). 

[Formula 32] 
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iY0RS=((MAX_KAICHO+1)/2*1000 + 299*iR0RS+587*iG0RS+1 14*iB0RS 
+ 500)/(1000) - (MAX.KAICHO+1 )/2 
= ( 

( 

((MAX_KAICHO+1)/2 « bitSHIFT) 
+ ((299«bitSHIFTH50O)/1O0O*iR0RS 
+ ((587<<bitSHIFT)+500)/1000*iGORS 
+ ((1 14«bitSHIFT)+500)/1000*iBORS 
+ (1 « (bitSHIFT-1)) 
) » bitSHIFT 
) -(MAX.KAICHO+D/2; 
iCbRS=((MAX_KAICHO+1)/2*2*886-299*iR0RS-587*iG0RS+(886) 
*(iB0RS+1))/(2*8B6)-(MAXJ<AICHO+1)/2 
= ( 

( 

((MAX.KAICHO+1 )/2 « bitSHIFT) 

- ((299«bitSHIF1>886)/(2*886)*iR0RS 

- «587«bitSHIFD+886)/(2*886)*iG0RS 
+' ((886«bitSHIFD+886)/(2*886)*iB0RS 
+ (1 « (brtSHIFT-1)) 

) » bitSHIFT 
) - (MAX_KAlCHO+1)/2; 
iCrRS=((MAX.KAICHO+1)/2*2*701~114*iB0RS-587*iG0RS+(701) 
*(iR0RS+D)/(2*701) - (MAX_KAICHO+l)/2 

' ( 

((MAX.KAICHO+1 )/2 « bitSHIFT) 

+ ((701<<b»tSHlFT)+701)/(2*701)*iRORS 

- ((587«bitSHIFT)+701 )/(2*701 )*i GORS 

- ((587«bitSHIFT)+701 )/(2*701 )*iGORS 
+ (1 « (bitSHIFT-1)) 

) » bitSHIFT 
) - (MAX_KAICHO+1)/2; 

The formula 32 can be simply expressed as 

follows : 
[Formula 33] 



ixORS= ( (1225*iR0RS + 2404*iGORS + 467*iB0RS 
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+ (257 « 11)) » 12) - 128; 
iCbRS= ((-691*iR0RS - 1357*iG0RS 

+ ((iBORS + 257) « 11)) » 12) - 128; 
iCrRS= ( ( ( ( iRORS + 257) « 11) 

- 1715*iG0RS - 333*iB0RS) » 12) - 128; 

The formula 33 can be further simplified as 

follows ; 
[Formula 34] 

iYORS= ( ( (129«12)+1225«iR0RS+2404KiG0RS+467xiB0RS 

+ <1«11))>>12)-128; 
iCbRS= ( ( <128«12) -691xiRORS-1357xiG0RS+2O4BxiBORS 

+ <1«11) )»12) -12B; 
iCrRS= ( ( (12B«12) -204BxiR0R5-1715*iG0RS-333xiB0RS 
+ <1«11)>»12) -128 ; 

In addition/ the formula 34 can be expressed 
by using powers of 2 as follows : 
[Formula 35] 



Y ( < (l28 x2 12 )+1225 xR(o)+2404 xG(o)+467 xB(o)+2 n ng 

Cb (o ) = ( 12S x 2 ' 2 )- 691xR ( 0 )- 13 ^ 7 xG (°K 2 " xB(o)+2" _ 128 . 

2 

Cr(o)^ ^ 12S x2l2 )~ 2U xR ^°^ 17 ^ 5 xG (°)- 333 xB(o)+2 u 12g 



2 l 

Regarding this color conversion, an available 
range of the input R(o)G(o)B(o) data is limited to realize 
the reversibility and the compatibility. That is , the 
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available range of the input R(o)G(o)B(o) data inputted in 
the formula 29 is limited as follows: 
[Formula 36] 

1SG(o) , and 

l^B(o) 

It is required to change values being out of boundary values 
to be the boundary values. If a value being out of the 
boundary value is input and the backward transformation is 
conducted in accordance with the formula 38 , a result from 
the formula 29 described above may be out of the available 
range (limitation range) . In this case, a problem occurs 
in that Y(o) , Cb(o) , and Cr(o) data can not be restored. 

In a case of restoring original Y (o) Gb (o) Cr (o) 
data, as the available range (limitation ragne) of the input 
R(o)G(o)B(o) data, the R(o) G (o)B (o) data being already 
obtained are used for formula 38 described later. And in 
a case of simply restoring Y (o) Cb (o) Cr (o) data, there is 
no special limitation , and restored Y (o) Cb (o) Cr (o) data can 
be obtained by using the R(o)6(o)B(o) data obtained from 
the following formula 37 describe later as the input data 
for the above converting formula. 

b. compatible conversion from integer values of the 
Y(o)Gb(o)Cr (o) data to integer values of the R(o)G(o)B(o) 
data 
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The backward transforming can be conducted by 
an integer operation of a color converting function shown 
by the following formula 37: 
[Formula 37] 

R(o) = 



0(o). 



2x(DxY(o) + 2x(P-x M )xCr(o))+P 
2xD 

px((D-x M -yJxDxY(o)-2xy M x(D-y M )xCb(o)-2xx M x(D^x M )xCr(o))>j 

[ + (P-x„-y M )xD J 

2x(D-x M -y M )xD 



b(o)=\ 2x ( DxY (°) +2x ( p -yM) x c b (°)) +p 

* ' L 2xD 

For example , similar to the formula 31/ in 

accordance with the JPEG 2000 algorithm/ when 

x=*299/1000=weighted coefficient of the red color , 

y=114/1000=weighted coefficient of the blue color, 

l-x-y=587/1000=weighted coefficient of the green color, 

and MAXkgb=256 levels, the formula 61 concerning the 

Y (o)Cb (o) Cf (o) data -> R(o) 6 (o)B (o) data conversion can be 

replaced with a color converting function shown by the 

following formula 38: 

[Formula 38] 

, x | (1000xY(o)+1402xCr(o))+500 I 

R W"L iooo J 

{ , I (587 xlOOO x V(o)-2 x 1 14 x886xCb(o)- 2 x299x701xCr(o))+ 587x500 I 

587x1000 J 
( x I (l 000 xY(o)+ 1772 xCb(o))+ 500 I 
W "L 1000 J 

In addition/ by applying a bit shift to the 
color converting function shown by the formula 38 , it is 
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possible to improve an operation speed. For example , since 
the level number of the RGB is 256 and an available range 
is from 0 to 255, the formula 31 can be expressed in C 
language where MAX_KAICHO=155 as a m a x i mum level value, 
bitSHIFT=12 (the number of bit shifts) as exponents of 
powers of 2, iRORS, iGORS, and iBORS denote R(o) , G(o) , and 
B(o) , and iYORS, iCbRS, and iCrRS denote Y(o) , Cb(o) , and 
Cr (o) . 
[Formula 39] 
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iR0RS=((MAX_KAICHO+1)/2*1000 + 1000*iYOS+2*701*»CrS + 500)/(1000) 

- (MAX KAICHO+D/2 

= ( 

( 

((MAXKAICHO+O/2 « bitSHlFT) 
+ ((1000«bitSHIFT)+500)/1000*iYOS 
+ ((2*701 «bitSHIFT)+500)/1000*iCrS 
+ (1 « (bitSHIFT-1)) 
)» bitSHIFT 
)-(MAX_KAICHO+1)/2; 

iG0RS=( 

(MAX .KAICHO+1 )/2*(2*5B7*1 000) 

+ 2*(587*1 000*iY0S-2*1 1 4*886*iCbS-2*299*701 *iCrS) 
+ 587*1000 
)/ (2*587*1000) 

- (MAX.KAICHO+D/2 
=( 

((MAX_KAICHO+1)/2 « bitSHIFT) 
+ ((2*587*1 000/4«bitSHIFT)+56 7*1 000/4)/(2*587*1000/4)*iY0S 

- ((114*8B6<<bltSHIFT)+587*1000/4)/(2*587*1000/4)*iCbS 

- ((299*701 «bitSHIFT)+587*1 000/4)7(2*587*1 000/4)*iCrS 
+ (1 « (bitSHIFT-1)) 

) » bitSHIFT 
) - (MAX.KAICHO+D/2; 
iB0RS=((MAX.KAICHO+D/2*1000 + 1000*iYOS+2*886*iCbS + 500)/(1000) 

- (MAX_KAICHO+1)/2 
=( 

( 

((MAX.KAICHO+O/2 « bitSHIFT) 
+ ((1000«bitSHIFT)+500)/1000*iYOS 
+ ((2*886«bitSH!FT)+500)/1 000*iCbS 
+ (1 « (bitSHIFT-1)) 
) » bitSHIFT 
)- (MAX.KAICHO+O/2; 

The formula 39 can be simplified as follows : 

[Formula 40] 
iR0RS=( ( ( (257 + <iY0S«l) ) « 11) 

+ 5743*iCrS) » 12) - 128; 
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iG0RS=<<(<257 + (iYOS«l)) « 11) 

- 1410*iCbS - 2925*iCrS) » 12) - 128; 
iB0RS«((((257 + (iTOS«l)) « 11) 

+ 7258*iCbS) » 12) - 128; 
The formula 40 can be further simplified as 

follows : 
[Formula 41] 
iR0RS-( ( (128«12) + 4096xiY0S + 5743xiCrS 

+ (1«11))»12) - 128; 
iG0RS=( ( (128«12) + 4096xiYOS) - HlOxiCbS - 2925*iCrS 

+ (1«11) )»12) — 128; 
iB0RS=( ( (128«12) +4096xiY0S + 7258xiCbS-»- (1«11))»12) 

- 128; 

The formula 41 can be expressed by using powers 
of 2 as follows: 
[Formula 42] 

p(o) (l28x2 ia ) + 2' 2 xY(o) + 5743 xCr(o) + 2 11 ^ 

2 

v x (128 x 2 12 )+ 2 n x Y(o)- 1410 x Cb (o)- 2925 x Cr(o)+ 2 11 m 
B , v (128 x 2 12 )+ 2 12 x Y(o)+ 7258 x Cb (o)+ 2 U m 

In order to realize the reversibility and the 
compatibility by operations in accordance with the 
above-described converting formulae, it is required to 
limit an available range of the input Y (o) Cb (o) Cr (o) data. 
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First, it is assumed that the input Y (o) Cb (o) Cr (o) data axe 
the Y(o)Cb(o)Cr (o) data obtained by the formula 29 
described above. In order for Y(o) "Cb(o) "Cr (o) " data in 
which an editorial process such as an image modification, 
an electronic watermark, or a like is additionally provided 
to the Y(o)Cb(o)Cr(o) data to be effective, based on a 
maximum displacement difference with respect to an original 
value, a boundary value is changed to be within an allowance 
range in the R(o)G (0)6(0) data, which are original for the 
Y (o) Cb (o) Cr (o) data, beforehand. 

In the following, in order to prevent image data 
from being out of the boundary value by the image 
modification, a method for changing the boundary value to 
be within the allowance range before the image modification 
will be described. Each influence degree of AY, ACb, and 
ACr of Y(o), Cb(o), and Cr(o) with respect to the 
R(o)G(o)B(o) data is shown by each of the following ratios: 

AR: AY: ACr =1 : 1 : (2 (D-3fe) /D) 

AG: AY: ACr: ACb=l:l: (2xm (D-x*) / ( <D-*M~y M ) *D) ) : 

<2y M <D-y M ) / ( (D-x«-y M ) *D) ) 

AB:AY : ACb =1 : 1 : (2 (D-y M ) /D) 
For example, when D=1000, xm=299, and yM=114, 

AR: AY: ACr -±71:1:1.40, 

AG: AY: ACr: ACb±?l:l: 0.71:0. 34, 

AB:AY :ACb*?l:l :1.77. 
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For example, with respect to G data, a change in increments 
of AY=1, A Cb=l , and A Cr=l becomes a displacement of 
2. 05 (=1+0. 34+0. 71) . Thus, in a case of MAXrgb=256 levels, 
before the editorial process is started, an allowance range 
of the green color is changed to 4^G^251. Accordingly, 
even if the displacement is newly added, it is possible to 
suppress within an original allowance range 1^6^254. 

By using the R(o) G (o) B (o) data where the change 
concerning the boundary values is conducted as described 
above, Y(o) , Cb(o) , and Cr(o) are calculated in accordance 
with the converting formula shown by the formula 29. Even 
if the displacement within the allowance range is applied 
to each of the calculated Y(o) , Cb(o) , and Cr(o) and the 
RGB conversion is conducted in accordance with the 
converting formula shown by the formula 37, it is possible 
to restore the Y (o) Cb (o) Cr (o) data where the displacement 
of the allowance range is applied. 

These conversions are connected to be a 
Y<o)Cb(o)Cr (o) -> R(o)G(o)B(o) -> Y (o) ' Cb(o) ' Cr (o) ' data 
conversion. Therefore, it is possible to restore the 
Y (o) Cb (o) Cr (o) data. That is, 

Y(o) 'Cb(o) 'Cr (o) '=Y(o)Cb(o)Cr(o> can be realized. 

In detail, in a case of Xm=114, D=1000, and 
MAX RGS =256, the color conversion can be conducted as 
follows : 
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(1) allowance range of an original RGB to generate 
y <o)Cb (o) Cr (o) ; input of the converting formula shown by 
the formula 29 

l^R(o)^254, 1^6(0)^254, l^B(o)^254. 

(2) conversion result range of Y(o)Cb(o)Cr(o) 
generated in the allowance range shown in the above (1) : 
input of the converting formula shown by the formula 37 

l^Y(o>^254, -126^Cb(o) ^127, 
-126^Cr(o)^127, 

(3) conversion result range of the RGB generated by 
the converting formula shown by the formula 37 in the 
conversion result range shown in the above (2) : input of 
the converting formula shown by the formula 29 

lSR(o)^255, l^G(o)^255, l^B(o)^255, 

(4) conversion result range of the ¥ (o) Cb (o) Cr (o) 
restored by the converting formula shown by the formula 29 
in the conversion result range shown in the above (3) 

l^Y(o) ^254, -126^Cb<o)^127, 
-126^Cr(o> ^127. 

Table 3 shows color conversion results from the 
color converting function. 
[Table 3 J 
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example 1. 

example 2. 

example 3. 

example 4. 
example 5. 



YCbCr 
data 



1 (Y) 

O(Cb) 

O(Cr) 



128(Y) 
O(Cb) 
O(Cr) 



105(Y) 
84 (Cb) 
106(Cr) 



112(Y) 
11 (Cb) 
-5(Cr) 
254 (Y) 
O(Cb) 
0(Cr) 



R(o)G(o)B(o) 
data 



KG) 
KB) 
KR) 



128(G) 
128(B) 
128(R) 



0(G) 
254(B) 
254 (R) 



112(G) 
131 (B) 
105(R) 
254(G) 
254(B) 
254 (R) 



restored 
■> YCbCr 
data 

KY) 

0(Cb) 

0(Cr) 



128 (Y) 
0(Cb) 
0(Cr) 



105(Y) 
84 (Cb) 
106 (Cr) 



112 (Y) 
11 (Cb) 
-5(Cr) 
254 (Y) 
0(Cb) 
0(Cr) 



According to the color conversion results , the 
YCbCr data are restored, that is, Y=Y' , Cb-Cb' , and Cr=Cr' 
are maintained. Therefore, it can be seen that it is 
possible to guarantee the data compatibility 
(reversibility) concerning the color conversion. 

As a result, for example, even if a color 
conversion is generally defined as a non-reversible 
transforming process in the JPEG 2000 algorithm 
schematically shown in FZ6.2, it is possible to change the 
reversible color conversion, it is possible to change the 
color conversion to be reversible by applying the 
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embodiment of the present invention. FIG. 2 illustrates a 
compressing/extending system similar to FIG.l. in FIG, 2, 
in a compression system, after the RGB -> YCbCr conversion 
is conducted by the forward transformation o£ the color 
space having minimum error, the lossless compressing 
process such as a wavelet conversion, an entropy coding , 
or a like is conducted, and then encoded data are stored 
in a memory and an external storage device or distributed 
through a communication path such as the Internet. In an 
extension system, the encoded data are read from the memory 
or provided through the communication path, and the encoded 
data by being compressed are reversibly extended in a 
reverse procedure. Consequently, the YCbCr -> RGB 
conversion is conducted by the backward transformation of 
the color space having the minimum error is conducted so 
as to recover image data from the encoded data. 

The reversible transformation can be realized 
by applying the color conversion described above to the 
forward and backward transformation of the color space 
having the minimum error. Therefore, even in JPEG 2000 
algorithm, the color conversion can be the reversible 
transformation, and an image cjuality can be guaranteed. 
Also, an analog type compatible color converting algorithm 
and a digital type compatible color converting algorithm 
can be realized. 
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In the following , a usage requiring the 
reversible color conversion will be described with 
reference to FIG. 3, as an application example in a case in 
that an editorial process such as a manipulation or a like . 

(1) First , a procedure for generating 4:2:2 YCbCr data 
from RGB based data will be described. The procedure 
realizes a perfect restored YCbCr data conversion of 4:2:2 
YCbCr -> RGB -> 4:2:2 YCbCr conversion 

When ROGOBO andRlGlBl are assumed as two pixels 
being adjacent each other as the RGB based data, Cb and Cr 
are calculated by using the color converting function shown 
by the formula 18 based on R= (R0+R1) /2 , G=(G0+G1) /2, and 
B=(B0+Bl)/2. YO and Yl are calculated by the color 
converting function shown by the formula 18 by the color 
converting function shown by the formula 18 by using ROGOBO 
and R1G1B1, respectively. 

ROGOBO and R1G1B1 are recalculated by 
conducting the backward transformation based on YO , Yl , Cb/ 
and Cr calculated above . If each of the recalculated ROGOBO 
and R1G1B1 is not greater than or equal to 0 and less than 
or equal to 255 f YO and Yl are recalculated based on the 
above R, G, and B (each average of ROGOBO and R1G1B1) . 

(2) Next, input (arbitrary) 4:2:2 YCbCr based data 
-> RGB -> 4:2:2 perfect recovered YCbCr data conversion is 
realized. A procedure for changing input 4:2:2 YCbCr based 
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data will be described. 

R, 6, and B can be arbitrary values 
independently from each other within a range O^R, G, B 
S»255. However , in addition to applying limitations of 16 
^Y^235, 16^Cb^240, and 16^Cr^240 to the YCbCr data, 
another condition is applied in that when R, G, and B are 
calculated from Y, Gb , and Cr, it is not effective if the 
conversion results of R, G, and B are not within ranges 0 
^R, G, B^255. Therefore, first , by using YCbCr - -> RGB 
converting expression by the color converting function 
shown by the formula 22 to the input 4:2:2 YCbCr based data, 
ROGOBO and R1G1B1 of the RGB based data described above (1) 
are calculated, in this case, when each calculation result 
becomes out of the boundary value (0,255) , the calculation 
result is changed to be the boundary value . 

Therefore f since ROGOBO and R1G1B1 become 
values being in the boundary values, the 4:2:2 YCbCr, in 
which ROGOBO and R1G1B1 are applied to a method described 
the above (1) and calculated, becomes the perfect 
reversible YCbCr. A stable value for the color conversion 
described in FIG. 3 is a value having reversibility in that 
an original value is restored when the backward 
transformation is further conducted to a value where the 
color conversion is conducted. 

(3) Next, a method for guaranteeing 4:2:2 YCbCR based 
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data ->RGB -> 4:2:2 perfect restored YCbCr data conversion 
where the editorial process (image modification, 
electronic watermark, or a like) is conducted will be 
described* 

ROGOBO and R1G1B1 of the RGB baaed data 
described above (1) are calculated by using the YCbCr -> 
RGB conversion shown by the formula 18 for the 4:2:2 YCbCr 
based data before the editorial process is conducted. When 
each conversion result becomes out of the boundary value, 
the conversion result is changed to be the boundary value. 
In this case, an influence to an image value by the editorial 
process (image modification) is considered. Based on a 
maximum displacement difference with respect to a based 
value, the boundary values of ROGOBO and R1G1B1 are changed 
to be within an allowance range. 

A method for changing the boundary values into 
the allowance range before the image modification is 
conducted, in that even if the image modification is 
conducted, image data is not out of the boundary value. An 
influence degree of each of displacement values AY, ACb, 
and ACr of Y Cb Cr with respect to ROGOBO and R1G1B1 is 
expressed as following ratio: 

AR: A Y: ACr ^1:1.17:1.60 

AG: AY: ACr; Acb^lzl .17:0.82:0.40 

AB: A Y :ACb=l:1.17 :2.02 
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For example, a change in increments of Ay=1, ACb=l, and 
ACr=l for G data becomes a displacement of 2.39 
(-1.17+0.82+0.40) and influences. Accordingly, if an 
allowance for the G data is changed to 3^G^252 before the 
editorial process is started, the G data can be set not 
to be out of the original allowance range of 0^G^255 even 
if the displacement is newly added. 

By using the ROG0B0 and R1G1B1 where the 
allowance for the G data is changed, 4:2:2 YCbCr is 
calculated by the method described in the above (1) . Even 
if the color conversion is conducted by adding the 
displacement to the allowable range to the calculated 4:2:2 
YCbCr, the 4:2:2 YCbCR data where the displacement of the 
allowable range is added can be restored. 

[Reversible Color Conversion Using Bit Shifts] 

Next, a reversible color conversion (YCbCr -> 
RGB -> YCbCr) using bit shifts will be described. A process 
of an original RGB -> YCbCr -> RGB -> Y'Cb'Cr' -> R'G'B' 
conversion will be described in the following, 

TCbCr (iY0S,iCbS,iCrS) = Y'Cb'Cr' 
<iYORS,iCbRS,iCrRS) is realized, that is, iYOS^iYORS, 
iCbS-iCbRS, and iCrS-iCrRS are realized. In this case, 
RGB=R'G'B' is realized. 

The original RGB data within ranges O^R, G, B 
^255 are converted into ranges O^Y^ 2 55, -128^Cb^l28, 
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and -128 ^Cr ^128 by the YCbCR conversion, converted into 
ranges -l^R, 6, B^256 by the RGB conversion, and converted 
into ranges O^Y' ^255, -128^Cb' ^128 , and ~128^Cr'^128 
by the Y'Cb'Cr' conversion. 

On the other hand, as for a stipulated range 
of an effective color that can be presented by 8 bits/color, 

for RGB type color, the following ranges are defined: 
0^R^255, 0^G^255, and 0^B^255. 

and for YCbCr type color, the following ranges axe 
defined : 

0^Y^255, -128^Cb^+127, and -128 ^Cr^+127 . 

In these stipulated ranges, data ranges 
effective to reversibly convert are ranges l^R, G, B^ 
2S4 for the original ijnage RGB, ranges 1^Y^254, 
-126^Cb^l27, and -126^Cr^l27 for the YCbCR conversion, 
ranges O^R, G, B^255 for the RGB conversion, and ranges 
l^Y'^254, -126^Cb' ^127, and -126SCr' ^127 for the 
Y'Cb'Cr' conversion. 

That is, the original image RGB within the 
ranges l^R, G, B^254 can be reversibly converted. 

However, the original image RGB when R, G, or 
B-0 or R, G, B=255 , YCbCR and RGB may be out of the stipulated 
ranges, respectively in the original RGB -> YCbCr -> RGB 
-> Y'Cb'Cr' -> R'G'B' conversion. For example, there are 
one case of Cb=128 or one case of Cr=128, and many cases 
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of R, G, 8»(-l) or 256. 

The RGB -> YCbCr conversion is conducted after 
values of the original image RGB are changed to be in 1 to 
254, respectively, so as that conversion results can be 
within the stipulated range. For example, a conversion 
table for converting into a predetermined value may he 
prepared. 

Next, a memory capacity for maintaining the 
conversion table will be described by considering the 
boundary values during the converting process , with respect 
to the following state A, state B, and state C. 
- State A 

The boundaries (Y<0 or 254 ^Y, Cb^-127 or 
127<Cb,Cr^-127 or 127<Cr) of the original image RGB -> 
YCbCr conversion are concerned. Cases of invalid data 
being out of the stipulated range are shown as follows. 
[Formula 43] 

Y < 0:0 case/ (256 x 256 x 256) color- 0% 

Y > 255 : 0 case/ (256 x 256x256) color =0% 
Cb < -128 : 0 case/ (256 x 256X 256) color =0% 

Gb > 127 : 1 case/ (256 x 256x256) color =6 x10"«% 

Cr < -128 : 0 case/ (256 X 256 X 256) color =0% 

Cr > 127 .- 1 case/ (256x256 x 256) color =6x1 0"«% 

- State B 

Only the YCbCr -> RGB conversion is conducted 
(no compensating process is conducted) . Then, similar to 
the above state A, the boundary values of the RGB are 
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considered and cases of the invalid data being out of the 

stipulated range are considered as follows . 

R < 0:11,783 cases/ (256x256x256) colors^O. 07% 
R < 255:11,883 cases/ (256x256x256) colors^O . 07% 
6 < 0: 6,171 cases/ (256x256x256) colors^O. 04% 
6 < 255: 6,117 cases/ (256x256x256) colors^O . 04% 
B < 0:14,408 cases/ (256x256x256) colors^O . 09% 
B < 255:14,529 cases/ (256x256x256) colors^O . 09% 

Accordingly, total 64,890 cases/ (256x256x256) colors*? 

0.39% become the invalid data. 
- State C 

The compensating process is conducted when the 
YCbCr -> RGB conversion. Similar to the above the states 
A and B, the boundary values of the RGB are considered , and 
cases of the invalid data being out of the stipulated range 
and change methods for these cases are considered. Table 
4 shows results . 
[Table 4] 
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out of 
boundary value 


change 
method 


found 
when Y 
calculation 


found 
when Cb 
calculation 


found 
when Cr 
calculation 


total 


R< 0 


change to R= 0 


O(case) 


7(case) 


I.IIB(case) 


1,125 






l 


•7 
9 


I.IIO 


1,123 


G< 0 


G= 0 


70 


4 


83 


157 


G>255 


G=255 


67 


4 


81 


152 


B< 0 


B= 0 


10 


2.714 


10 


2,734 


B>255 


B=255 


10 


2,720 


10 


Z740 


total 




158 


5.456 


2,417 


B.031 



As shown in the table 4, in the state C, total 8, 031 cases 
become the invalid data. 

For data which can not be reversibly converted 
to be with in the stipulated range by the bit shifts in 
accordance with the formula 34 and the formula 40 , in a case 
of using the conversion table as an exception, instead of 
applying the conversion table to the state B (64 , 891 cases) , 
it is effective to apply the conversion table to the state 
C (8,031 cases) . Accordingly, it is possible to reduce the 
memory capacity for maintaining the conversion table in the 
state C to be one eighth the memory capacity in the state 
B. 

As described above, by applying the conversion 
table to the state C (8,031 cases), it is possible to 
realized the reversible conversion for 16,000,000 colors 
including the following examples 1 to 5 shown in table 3. 
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[Details] 

FIG • 4 schematically shows the conversion 
results of a case in that arbitrary values in 0 to 255 are 
input for the RGB based data when the analog type color 
conversion using the formula 16 and the formula 18 is applied 
<R=104, G=112, and B=134) . As a result, Y=Y' , Cb-Cb' , and 
Cr=Cr' are realized. 

In a case in that arbitrary values in 0 to 255 
are input for the RGB based data when the digital type color 
conversion is conducted in accordance with ITU-R BT.601 
using the formula 21 and the formula 23 (R^104, G=112, and 
B=134) , the conversion results are schematically shown in 
FIG. 5 . As a result, Y=*Y' , Cb=Cb' , Cr*=Cr' , R=R' , G=G' , B=B' 
are realised. 

In a case in that arbitrary values in 1 to 254 
are input for the RGB based data when the reversible color 
conversion is conducted in accordance with the CTPEG 2000 
using the formula 31 and the formula 37 (R=l , G=254, and 
B=254) , the conversion results are schematically shown in 
FIG .6 by comparing with a conventional not-reversible 
conversion as comparative information. While the 
convention not-reversible conversion is in a not-restored 
state, in a case of applying the present invention, Y=Y' , 
Cb=Cb' , and Cr=Cf , R'=R' , G'=G' , and B'=B' are restored. 
[Variation] 
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In the above explanation, only the color 
conversion is focused, In general , the present invention 
can be applied to a reversible data conversion in that the 
forward transformation and the backward transformation 
reciprocally convert data between unit systems having 
different resolution levels, a first unit system having a 
lower resolution revel is used as the common unit system 
when the forward transformation and the backward 
transformation are conducted, the integer operation is 
conducted between data in the first unit system having the 
lower resolution level and data in a second unit system 
having a higher resolution level being higher than that of 
the first unit system. 

In the following, it will be described that a 
perfect reversible conversion can be realized by a 
successive conversion of the first unit system having the 
lower resolution level -> the second unit system having the 
higher resolution level -> the first unit system having the 
lower resolution level. 

Two resolution levels a and 0 (where a>0>O. 
For example, in a case of mm <- -> BMCJ conversion, one inch 
is expressed by a unit system expressed by a =540:l/100mm 
unit and another unit system expressed by j3-12O0BMQ) are 
assumed. After data are converted by transforming a 0 
unit system to an a unit system, the converted data are 
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further converted by converting to be back to the $ unit 
system. In a case in that this conversion is connected by 
the following relationship: 
[Formula 44] 

■/J— >Q? converting function (argument:/? data) 
= L(2a x(£ data)+J_)/2/3 j 

1 

round up if 5 or greater and round down if less than 5 

■ CX-+0 converting function ( argument : a data) 
.= L(20 x(ot data)j±ja_)/2c* J 

1 

round up if 5 or greater and round down if less than 5 

the following relationship is obtained. 
a x ( B data ) / $ - 0 . 5 

including numbers after the decimal point 

: rational number 

< B .-> a converting function (argument : jS data) 
— - integer according to the above definition 

X (J3 data)/ B + 0.5 

including numbers after the decimal point 

: rational number 

< a x (B data) - 0.5 

including numbers after the decimal point 

: rational number 

< B -> a converting function (argument : 0 data + 1) 

— integer according to the above definition 
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x (J3 data + 1) / 0 + 0*5 

including numbers after the decimal point 

: rational number. 

Since the 0 -> a converting function 
(argument : 0 data + 1) - 0 -> a converting function 
(argument z 0 data) > a / 0 - 1 > 0 , even after data adjacent 
to original $ data (an integer incremented or decremented 
by 1 and an original integer) , there is not duplicated data 
and (0 data) - 1 

< {0 data) - 0.5 x 0/a - 0.5 

< & x $ -> a converting function (argument : 0 data ) 

/a - 0.5 

< a -> 0 converting function (argument : 0 -> a 
converting function (argument : 0 data) 

^ 0 x 0 ~> a converting function (argument : 0 data) 

/a + 0.5 

^ (0 data) + 0,5 z fl/a + 0.5 

< ( £ data) + 1 . 

That is, since a «-> 0 converting function 
showing an integer value (argument : 0 -> a converting 
function (argument : # data)) = ( 0 data) can be obtained, 
the original 0 data can be perfectly restored. 

As a result, in a case in that by using the second 
unit system of BMU as the common unit system, when between 
the first unit system of 1/lOOmm and the second unit system, 
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the data conversion of 

BMU data -> l/100mm data -> BMU data 
is conducted, the reversible conversion is conducted by the 
integer operation so that the BMU data are always identical . 
In order to realize this reversible conversion, an integer 
operation in accordance with the following formula 
[Formula 45] 

(1/1 OOnun data) = l(2540 x (BMU data) + 600)/l200 J 
(BMU data) = L(l 200 x(l/l 00mm data) + 1270)/2540 J 

is conducted* 

Calculations of powers of 2 in the formula 45 

can be changed to bit shifts, so that it is possible to 

improve a calculation speed. 

[Formula 47 J 

(1 /100mm data) = 

( 

L«2540<<bitSHlFT)+ 600)71200 I x (BMU data)+(1 «(bitSHrFT-1 )) 

)»bitSHIFT; 

(BMU data) = 

( 

L((1 200«bitSHlFT)+ 1 270V254D I X (1 /100mm data)+(1«(bitSHIFT-1)) 
)»bitSHIFT: 

In the formula 47, since underlined parts are 
fixed values, the underlined parts are calculated 
beforehand. Therefore, it is possible to calculate at 
higher speed. Moreover, in a case of the number of bit 
shifts (bitSHIFT) =12, it is possible to convert without 
any error in a range of 0^ (BMQ data) ^2400. 
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Accordingly, in particular , for a conversion 
from the BMU data using an inch unit (inch, yard, foot, or 
a like) and a meter unit using l/100znm data between a unit 
used in a personal computer and another unit used in a 
printer, by using the BMU data having the lower resolution 
as the common unit system, the 1/ 100mm data having the higher 
resolution can be sufficiently guaranteed. Therefore, it 
is possible to realize the reversible conversion for length 
data that can be perfectly restored. 

Color conversion results obtained by using the 
converting function are shown in table 5. 
[Table 5] 



BMU 

data 



1/1 00mm 
data 



restored 

BMU 

data 



example 1. 


0BMU 


0(1/1 00mm) 


0BMU 


example 2. 


1 


2 


1 


example 3. 


2 


4 


2 


example 4. 


5 


11 


5 


example 5. 


65 


138 


65 


example 6. 


66 


140 


66 


example 7. 


666 


1410 


666 


example 8. 


1199 


2538 


1199 


example 9. 


1200 


2540 


1200 



According to the color conversion results , th« 
BMU data is restored and the data compatibility 
(reversibility) concerning a length conversion can be 
guaranteed. 
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In the following, a color conversion (YCbCr -> 
RGB => YCbCr) method for further guaranteeing the 
reversibility will be described. 

[Color Conversion (YCbCr -> RGB => YCbCr) Method for 
Further Guaranteeing The Reversibility] 
(1) preparation 

(1-A) expression of a process for rounding down fractional 
part 

In general, a division process for dividing 
positive fixed point data by different positive fixed point 
data in a computer rounds down the fractional part. In 
mathematics, a result from rounding down the fractional 
part of a real number A is expressed by putting the real 
number A in Gaussian symbols. Accordingly , the division 
result of the fixed point data by the computer is shown by 
using a modulo (showing a remainder of the division and shown 
by "mod") as follows. 
[Formula 48] 

V I V— VmodU 
UJ U 

In the formula 48, a result from dividing fixed point data 
V by fixed point data U by the computer is expressed. 

(1-B) expression of a process for rounding up if 5 or 
greater and rounding down if less than S 

Similar to the above (1-A) , a result from 
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rounding up if 5 or greater and rounding down if less than 
5 is shown as follows, in a division result from dividing 
positive fixed data Y by positive fixed data x (Y/X) by the 
computer . 
[Formula 49] 

2Y + X I 2Y + X-(2Y+X)mod(2X) 



+ X I 2Y + X- 
iX J 



2X J 2X 
(1-C) expression of a process for rounding up if 6 or greater 
and rounding down if less than 6 

Moreover, in the formula 49, a result from 
rounding up if 6 or greater and rounding down if less than 
6 is shown as follows, in a division result from dividing 
positive fixed data Z by positive fixed data W (Z/W) by the 
computer . 
[Formula 50] 



2Z + W-1 
2W 



2Z + W-1 -(2Z+W— 1)Tnod(2W) 



2W 

Based on explanations in the above (1-A) , (1-B) , 
and (1-C) , a color converting method for further 
guaranteeing the reversibility will be described in the 
following. That is, a color conversion at one side is 
conducted by rounding up if 5 or greater and rounding down 
if less than 5, and another color conversion at another side 
(opposite direction of the color conversion) is conducted 
by rounding up if 6 or greater and rounding down if less 
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than 6, so that the reversibility can be perfectly 
guaranteed. 

(2) method for perfectly restoring data by a successive 
conversion of a lower resolution level -> a higher 
resolution level -> a lower resolution level 

Two resolution levels a and 0 are defined. A 
unit system is converted from a j9 unit system to an a unit 
system and the converted unit system is further converted 
by using data in the converted unit system to be back to 
the 0 unit system (where a>/J>0) . In this case of the mm 
O ->BMU conversion, one inch is expressed in the resolution 
level a^2540 : 1/100 mm unit or the resolution level 0 
n:1200BMU unit. 
[Formula 51] 

• /?■— ► a converting function ( argument : 0 data) 
= L(2Qf x(£ data)±JL)/2£ J 

i 

round up if 5 or greater and round down if less than 5 

_ 2 g(0 data)+ 0 - (2 a{0 data)+ 0 )mod(2 0) 

20 

•at—*/} converting function ( argument : a data) 
- L(2/8 x(a data) +g-1 )/2Qf J 

i 

round up rf 6 or greater and round down if less than 6 

2 0(a data)+ a - 1 - (2 0{a data)+ a — 1 )mod(2 a) 

2 of 

In a case in that a relationship shown in the above formula 
51 is satisfied, when 
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[Formula 52] 

a — » 0 converting function ( argument ; /3—>a converting function(argument : £ data)) 
_ 2 fl(g data) + a - 1 - (2 flq data) +g-1 )mod(2 g) 

' 2 ffl 2 a(0)+ j8 )mod(2 tf _., 

_ 2 orW+ 0 - (2 0*0)+ 0 ) mod(2 ^ ^ )mod(2ff) 

2£ 



2a 

<*(£)+ yS -(2 <*(/?)+ /8)mod(2 0)+ a - 1 ^ 
(2 cr(0) + 0 -<2 ff(0)+ fl)mod(2 #>+ a - 1 )mod(2 a) J 



(Za(0)+0 



Za 

-(.2a(.B)+0)mod(20)+a-l \ 
(0)+ 0 -<2 a(ff)+ fl)mod(2£)+ a -1 )mod(2 g)J 

= (/? cam) + £±«Zlr /-») Z <OH2 «-*» 

= (0 data) + F( or. 
the following formula 53 is obtained. 
[Formula 53] 

-i < zd£+.lz2) * F(a . , )s <£±lnl) < n 

Therefore, since a result of the above calculations is 
obtained by an integer , this method can be expressed by using 

F(a , J3) ♦ 

That is, since the a->fi converting function 
(argument : a -> 0 converting function (argument z J3 
data) ) = ( $ data) showing the integer value, it is possible 
to perfectly restore original /3 data and it is possible 
to guarantee the perfect restore even i£ <x = B . 
(3) color converting method for guaranteeing the 
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reversibility 

(3-A) color converting expression capable of re-display 
(precise equation by a rational expression) 

Based on the above [Application to Original 
color converting function Bing Based on Brightness and 
Color Difference : JPEG 2000] , an original color conversion 
can be expressed by a precise equation using rational 
numbers. 

A YCbCr -> RGB -> YCbCr converting expression 



is described as follow: 
[Formula 54] 



R" 




"Y" 




"V" 




"R" 


G 




Cb 




Cb* 


=A 


G 


B 




Cr 




Cr' 




B 



In the formula 54, a matrix is denoted by a letter A and 
is expressed as follows. 
[Formula 55] 

1 0 2(l-x) 

t _ 2y(1 -y) 2x(l-x) 

1-x-y 1-x-y 

1 2(l-y) 0 



x 1-x-y y 

x 1-x-y 1 

2(i-y) ~~2(l-y) 2 

1 1-x-y y 

2 ~2(l-xj 2(l^x) 



Accordingly, the reversibility can be guaranteed. 
(3-B) color converting expression by integers (fixed point 
calculation) 



(a) YCbCr -> RGB conversion where calculation of 
rounding up if 6 or greater and rounding down if less than 
6 

Based on the formula 54 and the formula 55 shown 
in the above (3-A) , by applying the formula 52 shown in the 
above (3-C) , RGB can be expressed by YCbCr as follows. 
[Formula 56] 



G = 



4(D-x)Cr+D- 1 - (4(P-x)CrH-D- 1 )rood(2D) 

2D 

Z-Zrnod(2D(D-x-y)) 



2D(D-x-y) 



_ _ 4(D-y)C b+D- 1 - (4(D-y)Cb+D-l )mod(2D) 

2D 

where Z=2D(D - x - y) Y-4(Dy-y 2 )cb-4(Dx-x 2 )Cr+D(D - x - y)- 1 

(b) YCbCr -> RGB conversion where calculation of 
rounding up if 5 or greater and rounding down if less 'than 
5 

Based on the formula 54 and the formula 55 shown 
in the above (3-A) , by applying the formula 49 shown in the 
above (3-B) , RGB can be expressed by YCbCr by substituting 
the formula 56 as follows* In the following expression, 
Xmod(x) = 0 to (x-1) , that is, a minimum value is 0 (zero) , 
a maximum value is (x-1) , and a mod calculation is expressed 
by "0 - (x-1)"- 
[Formula 57] 
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Y = 



2Rx+2G(D-x-yH-2By+D 
2D 



_ 2Rx +2G(D- x -y)+2By +D _ 0~(2D-1) 
~ 2D 2D 

f2D 2 -x-y-l 
= Y -I- -L- -x(0-(2D-l))-(0~(2D(D-x-y)-l)) 

^y(0-(2D-l))-D(0-(2D-l)) 

= Y + Z y (D, x, y) 
Wie above formula 57 is calculated. 
[Formula 58] 

- - 2D 2 -D _ _ v. 2D 2 -x-y-l 



From the above formula 58, since Y' is an integer, Zy becomes 
0 (zero) . Accordingly, Y=Y' can be realized. Similarly, 
[Formula 59] 

2(D-y)Max R0B -xR-(D-x-y)G+(P-y)B+(D-y) 



2(D-y) 

^2(D-y)Max ftC8 -xR -(D-x ~y)G +(D -y)B +(D-y ^ 
_[ -(0~(2(D-y)-l)) J 



—Max 



RGB 



Max. 



= Cb + 



2(D-y) rcb 
^4D(D-y)Max RGB +2D z -2Dy+x+y-D+1 > 
+x(0-(2D-l))+(0~(2D(D-x-y)- 1)) 
-(D-yX0-(2D-l))-2D(0-(2(D-y)-l))J 



4D(D-y) 

= Cb +Z„(D, x,y) 

the above formula 59 is calculated. 
[Formula 60] 

4D(D-y)-2D-x-1 , y 4D-1 

1 < 4D(D-y) * Z b ( D . x .y^-^5-< 1 



—Max 



ROB 
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From the above formula 60 , since Cb r is an integer, Zb becomes 
0 (zero). Accordingly, Cb=Cb' can be realized. Moreover, 
[Formula 61] 



Cr'= 



2(D-x)Max RCB +(D-x)R-(D-x-y)G-yB+(D-x) I Max 

2(D-x) J RGB 

2(D-x)Max R0B +(D-x)R-(D-x-y)Q\ 

. yB+ (D- x ) -(q~(2(D-x)-1))J _ 
~ 2(D=x) MaXRCB 



= Cr + 



'4D(D-x)Max RGB +2D 2 -2Dx+x+y-D+1 > 
-(D-xX0~(2D-l)>l-(0-(2D(D-x-y)-l)) 
+y(0-(2D- l))-2D(0~(2(D-x)~ l)) 



—Max 



RGB 



4D(D-x) 
= Cr +Z P (D, x, y) 
the above formula 61 is calculated. 
[Formula 62] 

< , ^Wp-i s Zr(D ,,, y)s ^L< 1 

From the above formula 62 , since Cr ' is an integer , Z r becomes 
0 (zero) . Accordingly, Cr=Cr' can be realized. 

In addition, the color conversion according to 
the present invention can guarantee the reversibility for 
image data having more than 256 levels. 

According to the present invention, in the 
forward transformation and the backward transformation, a 
unit system having a lower resolution revel is used as the 
common unit system and the data conversion is conducted by 
the integer operation between unit systems having different 
resolution levels, by using the common unit system. 



79 



Therefore, a unit system having a higher resolution level 
is sufficiently processed and then it is possible to realize 
the reversible conversion where data can be perfectly 
restored. 

In particular , regarding the digital color 
conversion using a quantisation of analog Y, Cb, Cr signals 
between the unit system for the color space RGB and the unit 
system for the color space YCbCr , the unit system for the 
color space YCbCr having a lower resolution level is used 
as the common unit system. Thus, since the unit system for 
the color space RGB is sufficiently processed/ it is 
possible to realize the reversible conversion in a color 
space where data can be perfectly restored. 

Moreover, it is possible to realize the 
reversible conversion in the color space where data can be 
perfectly restored by a predetermined color converting 
function. 

Furthermore, in a case in that a unit system 
for a print in such as a color space YMC is used, it is 
possible to realize the reversible conversion in the color 
space where data can be perfectly restored* 

In particular, regarding the color conversion 
in accordance with the international standard specifying 
a data format for converting an analog video signal into 
the digital data between the unit system in the color space 
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RGB and the unit system in the color space YCbCr, the color 
space YCbCr having a lower resolution level is used as the 
common unit system. Therefore, since the unit system in 
the color space RGB having a higher resolution level can 
be sufficiently processed, it is possible to realize the 
reversible conversion of the color space where data can be 
perfectly restored. 

In particular, regarding the color conversion 
of the original color based on the brightness and the color 
difference between the unit system in the color space RGB 
and the unit system in the color space YCbCr , the color space 
YCbCr having a lower resolution level is used as the common 
unit system- Therefore, since the color space RGB having 
a higher resolution level can be sufficiently processed, 
it is possible to realize the reversible conversion in the 
color space where data can be perfectly restored. 

Moreover/ in the available ranges for 
R(o)G(o)B(o) data and Y(o)Cb(o)Cr(o) data, a limited range 
is provided in that a maximum value and a minimum value are 
limited. Therefore, without being out of the range of the 
color converting function, it is possible to conduct the 
compatible conversion where data can be perfectly restored. 

In particular , regarding the length conversion 
between the BMU type data using an inch unit system and the 
1/100 mm type data using a meter unit system in the personal 
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comput r and the printer , the inch unit system for the BMU 
type data having a lower resolution level is used as the 
common unit system. Therefore, since the meter unit system 
for the 1/100 mm type data having a higher resolution level 
can be sufficiently processed, it is possible to realize 
the reversible conversion where length data can be 
perfectly restored. 

The present invention is not limited to the 
specifically disclosed embodiments , and variations and 
modifications may be made without departing from the scope 
of the present invention. 

The present application is based on the 
Japanese Priority Applications No .2003-013585 filed on 
January 22/ 2003 and No . 2004-13138 filed on January 21, 2004, 
the entire contents of which are hereby incorporated by 
reference . 



